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I. INTRODUCTION 



Entanglement or quantum correlations between many space-separated subsystems has been recognized as one of 
the most intrinsic properties of quantum mechanics and provides the basis for many genuine applications of quantum 
information theory. It is, then, quite natural to look for physical situations in which quantum entangled states are 
obtained. Most of the theoretical and experimental effort has so far been devoted to unveil physical realizations 
of quantum states describing two quantum correlated subsystems. The search for physical systems displaying clean 
three-party entanglement is not simple. In this paper, we shall analyze decays of particles as a natural scenario for 
fulfilling such a goal. More precisely, we shall show that the decay of ortho-positronium into three photons corresponds 
to a highly entangled state. Let us now review what entanglement can be used for and why it is interesting to look 
for quantum correlation between more than two particles. 

In 1935 Einstein, Podolsky and Rosen starting from three reasonable assumptions of locality, reality and 
completeness that every physical theory must satisfy, argued that quantum mechanics (QM) is an incomplete theory. 
They did not question quantum mechanics predictions but rather quantum mechanics interpretation (^]. Their 
argument was based on some inconsistencies between quantum mechanics and their local-realistic premises (LR) 
which appear for quantum states of bipartite systems, S Ti^ C3> T~td 2 - It was m 1964 when Bell || showed that 
any theory compatible with LR assumptions can not reproduce some of the statistical predictions of QM, using a 
gedankenexperiment proposed in Q with two quantum correlated spin-i particles in the singlet state 

| S >=^(|01>-|10». (1) 

In his derivation, as it is well-known, quantum correlations or entanglement have a crucial role. Actually, the singlet 
state is known to be the maximally entangled state between two particles. The conflict between LR and QM arises 
since the latter violates some experimentally verifiable inequalities, called Bell inequalities, that any theory according 
to the local-realistic assumptions ought to satisfy. It is then possible to design real experiments testing QM against 
LR (for a detailed discussion see [^|). Correlations of linear polarizations of pair of photons were measured in 1982 
showing strong agreement with quantum mechanichs predictions and violating Bell inequalities Q. Nowadays, Bell 
inequalities have been tested thoroughly in favor of QM [0 . 

More recently, it has been pointed out that some predictions for quantum systems having quantum correlations 
between more than two particles give a much stronger conflict between LR and QM than any entangled state of two 
particles. The maximally entangled state between three spin-^ particles, the so-called GHZ (Greenberger, Home and 
Zeilinger) state § 

|GHZ> = -^(|000> + |111», (2) 

shows some perfect correlations incompatible with any LR model (see [|J and also |j| for more details). It is then of 
obvious relevance to obtain these GHZ-like correlations. Producing experimentally a GHZ state has turned out to be 
a real challenge yet a controlled instance has been produced in a quantum optics experiment pOj. 
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Entanglement is then important for our basic understanding of quantum mechanics. Recent developments on quan- 
tum information have furthermore shown that it is also a powerful resource for quantum information applications. For 
instance, teleportation [jix|] uses entanglement in order to obtain surprising results which are impossible in a classical 
context. A lot of work has been performed trying to know how entanglement can be quantified and manipulated. Our 
aim in this paper consists on looking for GHZ-likc correlations, which are truly three-party pure state entanglement, 
in the decay of ortho-positronium to three photons. The choice of this physical system has been motivated mainly by 
several reasons. First, decay of particles seems a very natural source of entangled particles. Indeed, positronium decay 
to two photons was one of the physical systems proposed long time ago as a source of two entangled space-separated 
particles |l2] ]. On a different line of thought, some experiments for testing quantum mechanics have been recently 
proposed using correlated neutral kaons coming from the decay of a </>-meson Jl3| . In the case of positronium, three 
entangled photons are obtained in the final state, so it offers the opportunity of analyzing a quantum state showing 
three-party correlations similar to other experiments in quantum optics. 

The structure of the paper goes as follows. We first review the quantum states emerging in both para- and ortho- 
positronium decays. Then, we focus on their entanglement properties and proceed to a modern analysis of the 
three-photon decay state of ortho-positronium. Using techniques developed in the context of quantum information 
theory, we show that this state allows in principle for an experimental test of QM finer than the ones based on the 
use of the singlet state. We have tried to make the paper self-contained and easy to read for both particle physicists 
and quantum information physicists. The first ones can found a translation of some of the quantum information ideas 
to a well-known situation, that is, the positronium decay to photons, while the second ones can see an application 
of the very recent techniques obtained for three-party entangled states, which allow to design a QM vs LR test for a 
three-particle system in a situation different from the GHZ state. 



II. POSITRONIUM DECAYS 
A. Positronium properties 

Let us start reminding some basic facts about positronium. Positronium corresponds to a e + bound state. These 
two spin-i particles can form a state with total spin equal to zero, para-positronium (p-Ps), or equal to one, ortho- 
positronium (o-Ps). Depending on the value of its angular momentum, it can decay to an even or an odd number of 
photons as we shall see shortly. 

Positronium binding energy comes from the Coulomb attraction between the electron and the positron. In the 
non-relativistic limit, its wave function is ]l4| 



( 27r )3/2 V> J ( 27r )3/2 7r(l + oV 

where a = -2-, i.e. twice the Bohr radius of atomic hydrogen, and m is the electron mass. Note that the wave 

ma 1 J ° ' 

function takes significant values only for three- momenta such that p < ^ <C m, which is consistent with the fact that 
the system is essentially non-relativistic. 

The parity and charge conjugation operators are equal to 

Up = U C = (-l) i+S , (4) 

where L and S are the orbital and spin angular momentum. Positronium states are then classified according to these 
quantum numbers so that the ground states are 1 So, with J PC = h , for the p-Ps and 3 S\+ 3 Di, having J PC = 
for the o-Ps. 

Positronium is an unstable bound state that can decay to photons. Since a rt-photon state transforms as Uc\n"f} = 
(— l) n |ri'y) under charge conjugation, which is an exact discrete symmetry for any QED process such as the decay 
of positronium, we have that the ground state of p-Ps (o-Ps) decays to an even (odd) number of photons . The 
analysis of the decay of positronium to photons can be found in a standard QED textbook [ p"4| . Para-positronium 
lifetime is about 0.125 ns, while for the case of ortho-positronium the lifetime is equal to approximately 0.14 fis [ [l6| . 

The computation of positronium decays is greatly simplified due to the following argument. The scale which controls 
the structure of positronium is of the order of \p\ ~ am. On the other hand, the scale for postrinomium annihilation 
is of the order of m. Therefore, it is easy to prove that positronium decays are only sensitive to the value of the 
wave function at the origin. As a consequence, it is possible to factor out the value of the wave function from the 
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tree-level QED final state computation H] . A simple computation of Feymann diagrams will be enough to write the 
precise structure of momenta and polarizations which describe the positronium decays. Furthermore, only tree-level 
amplitudes need to be computed since higher corrections are suppressed by one power of a. Let us now proceed to 
analyze the decays of p-Ps and o-Ps in turn. 



B. Para-positronium decay 

Para-positronium ground state decays into two photons. Because of the argument mentioned above, the determi- 
nation of the two-photon state coming from the p-Ps decay is simply given by the lowest order Feynmann diagram of 
e + e _ — > 77. Since positronium is a non-relativistic particle to a very good approximation, the three-momenta of e + 
and e~ are taken equal to zero, and the corresponding spinors are replaced by a two-component spin. This implies 
that the tree-level calculation of the annihilation of p-Ps into two photons is equal to, up to constants, 

M (e+e- — > 77) ~ x+ f M 2 X- , (5) 

where (see fbifl for more details) x± 1S the two-component spinor describing the fermions, = X T i* J 2, and M2 gives 

M 2 = xe|)-fc7 2x2 = A(fci,Ai;fc2,A2)/2x2, (6) 

perm 

where e * = e *(ki,Xi) stands for the circular polarization vector associated to the outgoing photon i and ^2x2 
is the 2x2 identity matrix. More precisely, for a photon having the three-momentum vector k — \k\k — 
|fc|(sin(9cos(/>, sin # sin 0, cos 0), the polarization vectors can be chosen 

e (fc, A) = — ^= (cos#cos</> — iAsin</>, cos 9 sin <j> + zAcos0, — sm9) , (7) 
v2 

where A = ±1 and they obey 

fc-e(fc,A)=0 k x e (fc,A) = -iXe (k, A) 

e (ki,Xi) ■ e {k j: Xj) = -- (l - X l X J 'k i ■ k/j . (8) 

From the expressions of the polarizaton vectors and the three-momentum and energy conservation, it follows that the 
scalar term A is 

A(k,X 1 ;-k,X 2 ) = ~(X 1 + X 2 ), (9) 

and it verifies 

A(k,+l;-k,+l) = -A(k,-l;-k,-l) 

A(k, +1; —k, -1) = -A[k, +1; -k, -1) = 0. (10) 

The two fermions in the para-positonium ground state are in the singlet state, |S'=0, S Z =Q) = (|^, — |) — \ — \-, 5)), 
and then, using the previous relations for A and (||), the two-photon state resulting of the p-Ps desintegration is 

l^) = ^(l++)-|--)). (11) 

The two- photon state resulting from p-Ps decay is thus equivalent to a maximally entangled state of two spin-i 
particles. This is a well-known result and was, actually, one of the physical system first proposed as a source of 
particles having the quantum correlations needed to test QM vs LR O] . 
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C. Ortho-positronium decay 



The ground state of ortho-positronium has J PC — 1 and, due to the fact that charge conjugation is conserved, 
decays to three photons. Repeating the treatment performed for the p-Ps annihilation, the determination of the 
three-photon state resulting from the o-Ps decay requires the simple calculation of the tree-level Feynmann diagrams 
corresponding to e + e~ — ► 777. Its tree-level computation gives, up to constants, 

M(e+e- — > 777) ~ xt M zX-, (12) 

and the 2x2 matrix M3 is equal to H] 



M 3 = £ 

cyclic perm. 

where 



((e 2 • e 3 - & ' 4) e I + (e 2 • 4 + e 3 ' • (13) 



^ = fc 2 xe*. (14) 

Using (||) we can rewrite M3 in the following way 

M 3 = a-V{ki 1 \i;k 2l \2;h,h), (15) 

where 

? = ((Ai -A 2 )(A 2 + A 3 ) ?*(*!, Ai) (e *(fc 2 , A 2 ) • e**(fc 3 , A 3 )) 
+ (A 2 - A 3 )(A 3 + Ax) e*(k, A 2 ) (?*(k, A3) ■ ?*(k, 

+(A 3 -Ai)(Ai+A 2 ) e , *(fc3,A 3 )(e*(fci,Ai)-e'*(fe2,A2))) . (16) 

Notice that the helicity coefficient (Aj — Xj){Xj + A&) for the cyclic permutations oiijk explicitly enforces the vanishing 
of the (+ + +) and ( ) polarizations, 

V(k,+; k, +; k, +) = V{k u -; k, -; fcs, -) = 0. (17) 

On the other hand, the rest of structures are different from zero 

V(k,-; k, +; k, +) = 2 g* *(k u -)(i - fc 2 • fe) 

U(fc 1; +; k 2 , -, h, -) = 2 e*(k, +)(1 - fc 2 ■ fc 3 ), (18) 

and similar expressions for the other cyclic terms. 

The original e + in the ortho-positronium could be in any of the three triplet states. It can be shown, using 
( |l2| ) and (|l5|), that when the initial positronium state is \S—1, S z — 1) — | i, i), the decay amplitude is proportional 
to Vi + 1V2, while the same argument gives —Vi + 1V2 for |S=1, S z =— 1) = | — |, — |) and — V2V3 for |5=1, S z =0) = 
^(||, — |) + I — f> l))- Now, considering the explicit expressions of the polarization vectors (0), with 6 = \ without 



loss of generality, and (|18|), it is easy to see that the three-photon state coming from the o-Ps decay is, up to 
normalization, 

l^o(fci, kM) = (i - k ■ k)(\ + +-> + |- -+)) 
+(i-fc 1 .fc 3 )d + _+) + !-+_)) 

+(l-fc 2 .fc 3 )(| -++) + ! + —)), (19) 
when the third component of the ortho-positronium spin, S z , is equal to zero, and 

\Mk, k, k)) = (1 - k ■ k)(\ + +-) - | - -+)) 
+(i-k-k)(\ + -+)-{-+-)) 

+(i-k-k)(\ -++>-! + —», (20) 

when S z = ±1. 

The final state of the o-Ps decay is, thus, an entangled state of three photons, whose quantum correlations depend 
on the angles among the momenta of the outgoing three photons. For the rest of the paper we will consider the first 
family of states (S z — 0) although equivalent conclusions are valid for the second one. In the next sections we will 
analyze the entanglement properties of the states \ipo(ki, k 2 , ^3)), using some of the quantum information techniques 
and comparing them to the well-known cases of the singlet and GHZ state. 
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III. ENTANGLEMENT PROPERTIES 



The quantum correlations of the three-photon entangled state obtained from the o-Ps annihilation depend on the 
position of the photon detectors, i.e. on the photon directions we are going to measure. Our next aim will be to choose 
from the family of states given by ([l9|) the one that, in some sense, has the maximum amount of GHZ- like correlations. 
In order to do this, we first need to introduce some recent results on the study of three-party entanglement. 

The set of states \ipo(ki, A;2, ^3)) form a six-parameter dependent family in the Hilbert space H.2 ® Ti-2 ® Ti.2, so that 
each of its components is equivalent to a state describing three spin-i particles or three qubits (a qubit, or quantum 
bit, is the quantum version of the classical bit and corresponds to a spin-i particle). Two pure states belonging to 
a generic composite system 7i.® N , i.e. N parties each having a d-dimensional Hilbert space, are equivalent as far as 
their entanglement properties go when they can be transformed one into another by local unitary transformations. 
This argument gives a lower bound for the entanglement parameters a generic state \4>) S Hf N depends on. Since 
the number of real parameters for describing it is 2 + , and the action of an element of the group of local unitary 
transformations {7(2)® is equivalent to the action of U(l) x SU(2)® , which depends on 3N + 1 real parameters, 
the number of entanglement parameters is bounded by 2 N+1 — (37V + 1). For our case this counting of entanglement 
parameters gives six, since we have N = 3, and it can be proved that this is indeed the number of nonlocal parameters 
describing a state in H2 <8> H-2 <8 H2 [0 ■ 

The above arguments imply that six independent quantities invariant under the action of the group of local unitary 
transformations will be enough, up to some discrete symmetry, to describe the entanglement properties of any three- 
qubit pure state. Given a generic state \<f>) S Ttf 3 

\<t>) = Y. u ^ k ) U,k= 1,2, (21) 

where \i), \j), \k) are the elements of a basis in each subsystem, A, B and C, the application of three local unitary 
transformations U A , U B and U c transforms the coefficients tijk into 

t'ijk = ^ UiixUfpUtntaPT (22) 

From this expression it is not difficult to build polynomial combinations of the coefficient Ujk which are invariant under 
local unitary transformations [jlT| , p~8| . These quantities are good candidates for being an entanglement parameter. For 
example, one of these invariants is 

*^^U 1 j 1 k 1 ti 1 j 2 k 2 ti2j2k 2 ti 2 j 1 k 1 — ^ t (Pa)' (23) 

where pa = tvBc{\'P){4\) is the density matrix describing the local quantum state of A (and the same happens for 
B and C). In [fl8| the six linearly independent polynomial invariants of minor degree were found (a trivial one is the 
norm) and a slightly modified version of these quantities was also proposed in [|19| . In the rest of the paper we will 
not consider the norm, so the space of entanglement parameters of the normalized states belonging to H2 ® Ti-2 <8> H-2 
has dimension equal to five. 

A particularly relevant polynomial invariant is the so-called tangle, r, introduced in p0|| . There is strong evidence 
that somehow it is a measure of the amount of "GHZ-ness" of a state [fl9|-p2f . It corresponds to the modulus of the 
hyperdeterminant of the hypermatrix given by the coefficients Ujk pql , which from ( ^l|) corresponds to 

T {\4>)) — |Hdet(tjjfc)| = e iii2 e i3i4 e jlj2 e j3j4 e kik 3 £k 2 k4ti 1 j 1 k 1 ti 2 j 2 k2~ti3j3k 3 ti i j i k4 , (24) 

where eoo = en = and eoi = — eio = 1- This quantity can be shown to be symmetric under permutation of the 
indices i,j, k. 

Because of the interpretation of the tangle as a measure of the GHZ-like correlations, we will choose the position of 
the photon detectors, from the set of states (|l9|), the ones that are associated to a maximum tangle. In the figure [l| it 
is shown the variation of the tangle with the position of the detectors. It is not difficult to see that the state of ( Jl9| ) 
with maximum tangle corresponds to the case k\ ■ fc2 = k\ ■ k^ — k2 ■ k$ — — A, i.e. the most symmetric configuration, 
that we shall call "Mercedes-star" geometry. The normalized state obtained from ( ]I^ ) for this geometry is 

IV>> = ^= (I + +-> + ! --+> + | + -+) + | -+-) + | -++> + ! + — ». (25) 
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Note that the GHZ state has tangle equal to \, while the value of the tangle of (Eq) is lower, 



4 

r(m = ^- (26) 

It is arguable that a "Mercedes-star" geometry was naturally expected to produce a maximum tangle state. Indeed, 
GHZ-like quantum correlations do not singularize any particular qubit. 

Let us also mention that the state we have singled out has some nice properties from the point of view of group 
theory. It does correspond to the sum of two of the elements of the coupled basis resulting from the tensor product 
of three spin-i particles, \ <g> \ (g) \ , Q 

m = +h+\l-h), (27) 



y/2 V 2 2' '2' 2' 



where 

l l + 5 ) = ^ (l + + - ) - 
The quantum correlations of (E5h will be now analyzed. 



")) 

-))• (28) 



IV. USEFUL DECOMPOSITIONS 

In this section, the state (^5|) will be rewritten in some different forms that will help us to understand better its 
nonlocal properties. First, let us mention that for any generic three-qubit pure state and by performing change of 
local bases, it is possible to make zero at least three of the coefficients of (|2l]) fl9|]25|| . A simple counting of 
parameters shows that this is in fact the expected number of zeros. This means that by a right choice of the local 
bases, any state can be written with the minimum number of coefficients tyfc, i.e. we are left with all the non-local 
features of the state, having removed all the "superfluous" information due to local unitary tranformations. For the 
case of the state ( p5| ) it is easy to prove [^6| that it can be expressed as 

|V^> = ^= (|001) + |010> + |100)) + ^|111), (29) 

which is the minimum decomposition in terms of product states built from local bases (four of the coefficients tijk are 
made equal to zero). 

An alternative decomposition, that will prove to be fruitful for the rest of the paper, consists of writing the state 
as a sum of two product states. This decomposition is somewhat reminiscent of the form of the GHZ state, which is 
a sum of just two product states, and is only possible when the tangle is different from zero ||l9| , ^l| as it happens for 
our state (see Eq). The state then can be written as 



2 f (l\ { l\ fl 



a(|000) + \aaa)), (30) 



1* 



where |0) = land a = f ^ J . We omit the details for the explicit computation of this expression since they can 

be found in |l9|pl| ]. It is worth noticing that o-Ps decay is hereby identified to belonging to an interesting type of 
states already classified in quantum information theory pl[ . 

The above decomposition allows for an alternative interpretation of the initial state as an equally weighted sum 
of two symmetric product states. Note that the Bloch vector, n — (sin cos sin sin 0, cos 8), representing the first 
local spinor appearing in ( |30| ) is pointing to the z axis, i.e. h\ — (0, 0, 1), while the second is located in the XZ plane 
with an angle of 120° with the z axis, i.e. fi2 — (-^,0,— ^). By performing a new unitary transformation, ( |30"| ) can 
be written as 

i*> -I ((:)»(:)* CM:) • (:) » (•))• (31) 

where c = cos 15°, s = sin 15°. Now, the two Bloch vectors are in the XZ plane, pointing to the 8 = 30° and 8 = 150° 
directions. The GHZ state corresponds to the particular case c = 1 and s = 0. 
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V. QUANTUM MECHANICS VS LOCAL REALISM 



The quantum correlations present in some three-qubit pure states show, as it was mentioned in the introduction, 
a much stronger disagreement with the predictions of a local-realistic model than any two-qubit entangled state. In 
fact, contrary to the case of the singlet state, no LR model is able to reproduce all the perfect correlations predicted 
for the maximally entangled state of three qubits 0. The state emerging from o-Ps decay is not a GHZ state, 
although it has been chosen the one with the maximum tangle in order to maximize GHZ-like correlations. In this 
section we will show how to use it for testing quantum mechanics against local-realistic models, and then we will 
compare its performance against existing tests for the maximally entangled states of two and three spin-^ particles. 
We start reviewing some of the consequences derived from the arguments proposed in M. 



A. QM vs LR conflict 



Given a generic quantum state of a composite system shared by N parties, there should be an alternative LR 
theory which reproduces all its statistical predictions. In this LR model, a state denoted by A will be assigned to the 
system specifying all its elements of physical reality. In particular, the result of a measurement depending on a set 
of parameters {n} performed locally by one of the parties, say A, will be specified by a function a\({n}). The same 
will happen for each of the space-separated parties and, since there is no causal influence among them, the result 
measured on A can not modify the measurement on B. For example, if the measurement is of the Stern-Gerlach type, 
the parameters labeling the measurement are given by a normalized vector n and a\{n) = a are the LR functions 
describing the outcome. 

The LR model can be very general provided that some conditions must be satisfied. Consider a generic pure 
state belonging to H2 <8> Ti-2 <8> Ti-2 shared by three observers A, B and C, which are able to perform Stern-Gerlach 
measurements in any direction. Since the outcomes of a Stern-Gerlach measurement are only ±1, it is easy to check 
that for any pair of measurements on each subsystem, described by the LR functions a and a' , b and b' , c and c', and 
for all their possible values, it is always verified 

a'bc + ab'c + abc' - a'b'c' = ±2. (32) 

It follows from this relation that 



- 2 < (a'bc + ab'c + abc - a'b'c) < 2. (33) 

This constraint is known as Mermin inequality p7[ and has to be satisfied by any LR model describing three space- 
separated systems. 

Let us now take the GHZ state (|^). It is quite simple to see that if the observables a and a' are equal to a y and 
a x (the same for parties B and C), the value of ( |33"| ) is —4, so an experimental condition is found that allows to 
test quantum mechanics against local realism. Note that this is the maximal violation of inequality ( |3"5| ) . Moreover, 
the GHZ state also satisfies that a'bc — ab'c = abc' = —a'b'c' = — 1 and no LR model is able to take into account 
this perfect correlation result because of (|3^) ||. This is a new feature that does not appear for the case of a two 
maximally entangled state of two spin-^ particles. In this sense it is often said that a most dramatic contrast between 
QM and LR emerges for entanglement between three subsystems. 

Let us go back to the state given by the ortho-positronium decay (p5|). Our aim is to design an experimental 
situation where a conflict between QM and LR appears, so we will look for the observables that give a maximal 
violation of (|33|). Such observables will extremize that expression. Using the decomposition (|3l|), the expectation 
value of three local observables is 



(abc) = (ip\(h a ■ a) ® (n b ■ a) ® (n c ■ a)\ip) 

\ i— a,6,c i— a.fc.c 



Y\_ (c cos 6*j + s sin 6*j cos^j) + J~J (— ccos#j + ssin#j cos^>j) 



Yl sin^(c 2 e-^ +s 2 e^)+ JJ sin0 4 (cV^ + s 2 ^) ] , (34) 

i— a,b,c i—a,b,c 
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where c = c 2 — s 2 and s = 2sc. Because of the symmetry of the state under permutation of parties, the Stern- 
Gerlach directions are taken satisfying h a = hb — h c = (sin 9 cos 4>, sin 9 sin 0, cos 9) and h a > = hy = h c i = 
(sin 9' cos <fi', sm9' sin</>', cos 9'). Substituting this expression in (5^), we get the explicit function f(9,(p,9\4>') to 
be extremized. For the case of the GHZ state described above, the extreme values were obtained using two observ- 
ables with 6 = 6' = ^, i.e. in the XY plane. Since ( |3l| ) is the GHZ-like decomposition of the initial state, we take 
9 = 9' = ^ and it is easy to check that in this case ^\ g _ g ,_„ = §gr\g- g i-iL — 0, V<f>,<fi'. Mantaining the parallelism 

with the GHZ case, it can be seen that all the partial derivatives vanish when it is also imposed <i> — \ and 4>' = 0. In 
our case the calculation of (^) gives —3, so a conflict between local-realistic models and quantum mechanics again 
appears, and then the three-photon state coming from the ortho-positronium decay can be used, in principle, to test 
QM vs LR with the set of observables given by the normalized vectors 

h a = fib = h c = (0, 1,0) h a ' = fib' = h c ' = (1,0,0). (35) 

There is an alternative set of angles <fr and <j)' that makes zero all the partial derivatives of /: the combination of local 
observables ( |3^ ) is equal to w —3.046 for 

VI7+I^ 24 . (36) 

This second set of parameters will be seen to produce in the end a weaker dismissal of LR. 

Our next step will be to carry over the comparison of this QM vs LR test against the existent ones for the maximally 
entangled states of three and two spin-i particles, i.e. the GHZ and singlet state. It is quite evident that the described 
test should be worse than the obtained for the GHZ state. It is less obvious how this new situation will compare with 
the singlet case. 



B. Comparison with the maximally entangled states of two and three spin-i particles 

We will now estimate the "strength" of the QM vs LR test proposed above, being this "strength" measured by 
the number of trials needed to rule out local-realism at a given confidence level, as Peres did in |p8[. A reasoning 



anologous to the one given in |28j will be done here for the state (25) and the observables (p5| 



Imagine a local-realistic physicist who does not believe in quantum mechanics. He assigns prior subjective proba- 
bilities to the validity of LR and QM, p r and p q , expressing his personal belief. Take for instance = 100. His LR 
theory is not able to reproduce exactly all the QM statistical results of some quantum states. Consider the expectation 
value of some observable O with two outcomes ±1 such that (O) = E q is predicted for some quantum state, while 
LR gives (O) — E r ^ E q . Since the value of the two possible outcomes are ±1, the probablity of having O = +1 is 
q = 1+ 2 q for QM and r = 1+ ^ r for LR. An experimental test of the observable O now is performed n times yielding 
m times the result +1. The prior probabilities p q and p r are modified according to the Bayes theorem and their ratio 
has changed to 

Pr _ Pr p{m\ L R) , 3? s 



Pq PqP{rn\ QM ) 
where 



id 



p(m\ LR )=(Z)r m (l-rr-™, (38) 



is the LR probability of having m times the outcome +1, and we have the same for p{tti\qm)i being r replaced by q. 
Following Peres p8| , the confidence depressing factor is defined 

Ds *M«») (ly(lzl)™ , (39) 



which accounts for the change in the ratio of the probabilities of the two theories, i.e. it reflects how the LR belief 
changes with the experimental results. Like in a game, our aim is to destroy as fast as we can the LR faith of our 
friend by choosing an adequate experimental situation. It can be said, for example, that he will give up when, for 
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example, D — 10 4 . Since the world is quantum, m = qn, and the number of experimental tests needed to obtain 
D = 10 4 is equal to 



n D(q, r) = — j- r- = — -, (40) 



being K(q,r) the information distance |29] between the QM and LR binomial distribution for the outcome +1. The 
more separate the two probability distributions are, measured in terms of the information distance, the fewer the 
number of experiments njj is. 

Let us come back to the three-party entangled state coming from the ortho-positronium decay ( |25| ) under the local 
measurements described by (|35|). As it has been shown above, a contradiction with any LR model appears for the 
combination of the observables given by the Mermin inequality In our case quantum mechanics gives the following 
predictions 

(a'bc) = (ab'c) = (abc') = -\ (a'b'c') = +1, (41) 

and this implies that q\ = prob(a'bc = +1) = prob(ab'c = +1) = prob(abc' = +1) = i and q 2 = prob(a'b'c' = +1) = 1. 
This is the QM data that our LR friend has to reproduce as well as possible. Because of the symmetry of the state he 
will assign the same probability r\ to the events a'bc — +1, ab'c — +1 and abc' = +1 and r 2 to a'b'c' = +1. However, 
his model has to satisfy the constraint given by (|33|), so the best he can do is to saturate the bound and then 

3n = r 2 => < n < i (42) 

Now, according to the probabilities r\ and r 2 his LR model predicts, we choose the experimental test that minimizes 
(fhj), i.e. we consider the event a'bc = +1 {a'b'c' = +1) when ni){qi,ri) < nr>(q 2 ,r 2 ) (^d(<Z1) ri) > n£,(q 2 ,r 2 )), and 
the experimental results will destroy his LR belief after nr>(qi,ri) (tid(<72> ^2)) trials. The best value our LR friend 
can assign to r\ is the solution to 

nD(qi,n) =n D (q 2 ,r 2 ), (43) 

with the constraint (^|), and this condition means that r x ps 0.315 and njj » 161 trials are needed to have a depressing 
factor equal to 10 4 . Repeating the same calculation for the observables giving by (|36|), the number of trials slightly 
increases, nu ~ 166, despite of the fact that the violation of the inequality is greater than the obtained for (|35|) . 

In ref. |28| | the same reasoning was applied to the maximally entangled state of two and three spin-i particles, 
showing that Ud ~ 200 in the first case, and ud ~ 32 for the latter (see table ||). Our result then implies that the 
three-photon entangled state produced in the ortho-positronium decay has, in some sense, more quantum correlations 
than any entangled state of two spin- ^ particles. 



C. Generalization of the results 



It is easy to generalize some of the results obtained for the entangled state resulting from the o-Ps decay. As it has 
been mentioned, this state can be understood as an equally weighted sum of two symmetric product states, since it 
can be written as (pl|). The Bloch vectors of the two local states appearing in this decomposition form an angle of 
120°. It is clear that the conclusions seen above depend on the angle between these vectors, i.e. with their degree of 
non-orthogonality. The family of states to be analyzed can be parametrized in the following way 

'*"» =-((::M:M:) + C:M:) «(:))■ w 

where S is the angle between the two local Bloch vectors, eg = cos (^j^) and sg = sin ( ZL j^) and as is a positive 
number given by the normalization of the state. An alternative parametrization of this family is, using ( p9| ) and 
defining 5' = |, 

\tf)(5)) = 2a s (sin 2 S'cosS' (|001) + |010) + 1 100)) + cos 3 <5'|H1)) • (45) 
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The expectation value of three local observables for this set of states follows trivially from ( |34| ) . Using this expression 
it is easy to see that the combination of the expectation values of ([33]) has all the partial derivatives equal to zero for 
the set of observables given in (|5|) independently of S. For these observables, the dependence of expression ([33]) with 
the degree of orthogonality between the two product states is given in figure |^. There is no violation of the Mermin 
inequality for the case in which 6 < 85°. In this situation one can always found a LR model able to reproduce the 



QM statistical prediction given by (|33j) and the observables (|35|). We can now repeat all the steps made in order 
to determine the number of trials needed to rule out local realism as a function of the angle S. In figure [| we have 
summaryzed the results. We have shown only the cases where the number of trials is minor than two hundred, since 
this is the value obtained for the singlet. Note that the case S — 120°, which corresponds to (p5f), is very close to the 
region where there is no improvement compared to the maximally entangled state of two qubits. 

All these results can be understood in the following way: the smaller the angle between the two local states, 5, the 
higher the overlap of the state with the product state having each local Bloch vector pointing in the direction 

of the x axis, which corresponds to the state |111) in (fl5|). This means that the quantum state we are handling is too 
close to a product state P5|, and thus, no violation of the Mermin inequality can be observed. 



VI. CONCLUDING REMARKS 



In this work we have analyzed the three-particle quantum correlations of a physical system given by the decay of the 
ortho-positronium into a three-photon pure state. After obtaining the state describing the polarization of the three 
photons (p5|), some of the recent techniques developed for the study of three-party entanglement have been applied. 
The particular case where the three photons emerge in a symmetric, Mercedes-star-like configuration, corresponds to 
the state with the maximum tangle. We have shown that this state allows a priori for a QM vs LR test which is 
stronger than any of the existing ones that use the singlet state. In this sense, ortho-positronium decays into a state 
which carries stronger quantum correlations than any entangled state of two spin-^ particles. 

Bose symmetrization has played a somewhat negative role in reducing the amount GHZ-ness of the o-Ps decay state. 

Indeed, the natural GHZ combination |+H — ) + | h) emerging from the computation of Fcynmann diagrams has been 

symmetrized due to the absence of photon tagging to our state |++— ) + H h) + | — ++) + | h) + |— +— ) + |H ), 

inducing a loss of tangle. The quantum optics realization of the GHZ state does avoid symmetrization through a 
geometric tagging pj. It is, thus, reasonable to look for pure GHZ states in decays to distinct particles, so that 
tagging would be carried by other quantum numbers, as e.g. charge. It is, on the other hand, peculiar to note that 
symmetrization in the K°K° system is responsible for its entanglement (| H — ) + | — h)) El]. 

Finally, let us briefly discuss the experimental requirements needed for testing quantum mechanics as it has been 
described in this paper. In order to do this, the circular polarizations of the three photons resulting from an ortho- 
positronium decay have to be measured. The positions of the three detectors are given by the "Mercedes-star" 
geometry and their clicks have to detect the coincidence of the three photons. The energy of these photons is of 
the order of 1 Mev. Polarization analyzers with a good efficiency would allow us to acquire statistical data showing 
quantum correlations which would violate the Mermin inequality discussed above. Unfortunately, as far as we know, 
no such analizers exist for this range of energies. A possible way-out might be to use Compton scattering to measure 
the photon polarizations fl30|| . However, Compton effect just gives a statistical pattern depending on the photon and 
electron polarizations which is not a direct measurement of the polarizations. Further work is needed to modify our 
analysis of QM vs LR to accommodate for such indirect measurements. 
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FIG. 1. Variation of the tangle with the position of the photon detectors, that are represented by two angles (the third one 
has to sum up to 360°). We have taken r = when the position of the detectors, i.e. the photon trajectories, are incompatible 
with momentum conservation. 



Violation of inequality 




Angle 



FIG. 2. Violation of the Mermin inequality ( j33[) with the angle 5 for the family of states (fl4|). We have substracted 2 to 
the combination of the expected values of (^3|), so a positive value means that a conflict between QM and LR appears. 
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FIG. 3. Number of trials needed to rule out local realism as a function of the angle 5 for the family of states ([44|). Values 
greater than two hundred are not shown since in these cases there always exists a two-qubit entangled state which gives the 
same result, i.e. it has the same "strength" for ruling out local realism. 
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State 


Number of trials 


GHZ 

Positronium state (Ba) 
Singlet 


fa 32 
« 161 
« 200 



TABLE I. Comparison of the strength of the QM vs LR test which can be performed for the maximally entangled states of 
two and three spin-| particles and for the three-photon entangled state reulting from the ortho-positronium annihilation. 
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